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Abstract 

We study the three-quark static potential in perturbation theory in QCD. A complete next-to- 
leading order calculation is performed in the singlet, octets and decuplet channels and the potential 
exponentiation is demonstrated. The mixing of the octet representations is calculated. At next-to- 
next-to-leading order, the subset of diagrams producing three-body forces is identified in Coulomb 
gauge and its contribution to the potential calculated. Combining it with the contribution of the 
two-body forces, which may be extracted from the quark-antiquark static potential, we obtain the 
complete next-to-next-to-leading order three-quark static potential in the colour-singlet channel. 

PACS numbers: 12.38.-t,12.38.Bx,14.20.-c,14.40.Pq 



1 



I. INTRODUCTION 



The interaction among heavy quarks has been explored since the QCD inception as an 
important tool to learn about the characteristics of the non-Abelian gauge dynamics in 
general and the QCD low-energy behaviour in particular jl, 3]. 

The QQ static potential is a very well known quantity for its crucial role in quarkonium 
phenomenology 0, 4 ] and for having been studied extensively by lattice gauge theories since 
their introduction 5). The typical shape of the colour-singlet QQ static potential, which is 
characterized by a short-range Coulomb behaviour and a long-range linear rise, well repre- 
sents the double nature of QCD as an asymptotically free and infrared confined theory Also 
gluonic excitations of a static quark-antiquark pair have been explored by lattice calcula- 
tions both in the long range, where they exhibit a stringy behaviour like the colour- singlet 
potential, and in the short ranee where they show a Coulomb-like behaviour in one of the 

nn 

two possible quark-antiquark colour configurations: singlet or octet [5|, |6[. 

More recently, non-relativistic effective field theories of QCD have provided a new way 
to look at the quark-antiquark potential and allowed, specially in the short range, calcu- 
lations with unprecedented precision 7j. Presently, the static quark-antiquark potential is 



complete* know, up to two loops W . Starting from three .oops the potential exhibits in- 
feared divergences; these have been ealcu,ated at leading order (LO) (itj and next-to-leading 
order (NLO) [111, a nd resummed at leading logarithmic (LL) [12| and next-to-leading log- 
arithmic (NLL) |l3] accuracy. The fermionic part of the three-loop finite contribution has 
been calculated recently [ijj]. High-order perturbative calculations show a remarkably good 
agreement with the lattice determinations of the static quark-antiquark energy up to a dis- 



15Ml7j. which allows to constrain the size of the unknown 



tance of about 0.2 - 0.3 fm Il3l 
higher-order contributions. 

The static quark-antiquark energy may be extracted from the large-time behaviour of the 
static quark-antiquark Wilson loop. Extremely accurate lattice determinations of the static 
energy at short distances (the smallest distance being about 0.08 fm) can be found in fl8 |. 
Also gluonic excitations between static quark-antiquark sources have been explored in the 
framework of effective field theories 19( and by means of lattice calculations [6|. Again, 



high-order perturbative calculations show agreement with accurate short-range lattice data 
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and allow for the precise extraction of the so-called gluelump masses 20] . Short distance 
studies of the quark-antiquark interaction tell us about the interplay of perturbative and non- 
perturbative contributions in QCD, in particular that perturbative contributions describe the 
data with a high accuracy up to distances of 0.2 - 0.3 fm, while a confining string sets in only 
at distances of about 0.5 fm 2l|], and that the operator product expansion does not appear 
to be violated. It is only natural to ask if these features are specific of quark-antiquark 
systems, i.e. mesons, or may also show up, and, in case, to which extent, in three-quark 
systems, i.e. baryons. 

The potential that describes the interaction of three heavy quarks Q is much less known 
than the heavy QQ potential, one of the reasons being the difficulty of producing QQQ states 
and the consequent lack of experimental data. This has led to a wide use of phenomenolo- 



gical models 
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23|, sometimes based on strong-coupling expansion arguments and lattice 



evaluations of the three-quark static Wilson loop; often a sum of two-body interactions has 
been used. 

A rigorous definition of the QQQ potential is provided by the non-relativistic effective 
field theory for QQQ states formulated in 24|] {QQq states have been considered in 24. l25|). 
This effective field theory is the three heavy-quark version of potential non-relativistic QCD 
(pNRQCD), the effective field theory first introduced for quarkonium in 19| |26j. pNRQCD 
is constructed from QCD as an expansion in the inverse of the heavy-quark mass m and 
in the distances between the heavy quarks (multipole expansion). At zeroth order in the 
multipole expansion, the equation of motion of pNRQCD is the Schrodinger equation with 
the potentials given by the Wilson coefficients of the six-fermion operators. The Wilson 
coefficients are calculated by equating, i.e. matching, amplitudes in QCD with amplitudes 
in pNRQCD order by order in 1/m and in the multipole expansion. In particular, the static 
potentials of the different colour representations are evaluated by matching to static Wilson 
loops in QCD. At distances shorter than the inverse of the typical hadronic scale, Aqcd, the 
degrees of freedom of pNRQCD are a QQQ colour-singlet field, two QQQ colour-octet fields, 
a QQQ colour-decuplet field, light quarks and low-energy gluons. The Wilson coefficients 
of the corresponding six-fermion operators are the singlet, octet and decuplet potentials 
respectively. They may be evaluated in perturbation theory. To the best of our knowledge 
only the LO expressions (excluding octet mixing) have been considered so far. At distances 
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larger than 1/Aqcd, when confinement sets in, the degrees of freedom of pNRQCD are only 
the QQQ colour- singlet field and light hadrons. Gluonic excitations of heavy-quark bound 
states cannot be resolved at such distances because of the mass gap of order Aqct) that they 
develop with respect to the colour-singlet state (cf. with the lattice data in [271, l28j). In 
this situation, the matching to pNRQCD cannot be performed in perturbation theory but 
must rely on non-perturbative methods. The non-perturbative static, spin-dependent and 
1/m colour-singlet QQQ potentials have been expressed in terms of Wilson loops in (for 



earlier work see 
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30] ). So far only the static potential has been evaluated on the lattice. 



Most of the existing lattice studies of the three-quark static potential have explored the 
region of large interquark distances 28|, I3ll438l]. As for the QQ case, the characteristic 
signature of the long-range non-Abelian dynamics is believed to be a linear "stringy" rising 
of the static interaction. Moreover, the general expectation for the baryonic case is that, 
at least classically, the strings meet at the so called Fermat (or Torricelli) point, which has 
minimum distance from the three sources (F-shape configuration). If this is the case, one 
should see a genuine three-body interaction among the static quarks. In another model 
39l | . the long range QQQ potential is simply the sum of two-body potentials (A-shape 
configuration). Most of the lattice calculations of the QQQ static potential have focused on 
distinguishing the Y configuration (favoured by data) from the A configuration, despite the 
difference between a A and a Y shape potential being rather small and difficult to detect. 
Recently, however, some data have accumulated that include short distances both at zero 
and finite temperature, and both for the lowest and for some higher gluonic excitations 
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40| . This opens the possibility to address, also for the QQQ system, questions 



about the short-range behaviour of the static potential and its gluonic excitations, and more 
specifically about the region of validity of perturbation theory and about the cross-over 
region from perturbative to non-perturbative QCD. In general, one expects this cross-over 
to happen in a more spectacular way than in the quark- ant iquark case, due to the overcoming 
of the long-range three-body forces over the short-range two-body Coulomb forces. 

In the paper, we focus on the potential between three static quarks in the different colour 
configurations and at short distances. Surprisingly, very little is known about it besides the 
LO expression. For all colour configurations, we will perform a complete NLO calculation 
showing explicitly how the exponentiation works at this order. For the singlet and decuplet 
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potentials, we will prove that the naive extension of the NLO two-body result turns out 
to be correct. For the octet potentials, we will need to account for the mixing, which 
already sets in at LO. At next-to-next-to leading order (NNLO) the first genuine three-body 
contribution appears. We calculate it for the singlet and decuplet colour configuration. In 
the colour-singlet case, combining the three-body contribution with the two-body one that 
can be extracted from the quark-antiquark static potential, we will obtain the complete 
NNLO potential. 

The plan of the paper is the following. In Sec. Ull we introduce the three-quark Wilson 
loop and define the potential. In Sec. Ill II we derive its expression at order g 2 for the singlet, 
the octets and the decuplet representations, showing that the two octet representations mix. 
In Sec. IIV1 we calculate the static potentials at order g 4 and show how exponentiation works 
at this order; a generalization of this result to N quarks in SU(iV) is provided in Sec. |Vj In 
Sec. IVTl we identify the first genuine three-body contribution to the potential that appears 
in perturbation theory at order g 6 and evaluate it in several geometrical configurations. In 
Sec. IVHl we derive the two-body colour-singlet contribution and hence provide the complete 
colour-singlet static potential at order g G . Sec. IVIIIl is devoted to the conclusions and a short 
outlook. Some technical details may be found in the appendices. 



II. THE THREE-QUARK STATIC POTENTIAL 

In this section, we consider the perturbative static potential of three heavy quarks. In the 
effective field theory language of 24|, the potentials in the different colour representations 
are the matching coefficients of the six-fermion operators made of two singlet, two octet 
or two decuplet fields. The matching coefficients can be ordered in powers of 1/m, the 
static potential corresponding to the first term in the series. The perturbative expression of 
the potential is expected to describe correctly the potential at short distances r, for which 
a s (l/r) < 1 holds. 

The static potential is computed by matching the appropriate Green's function in QCD 
with static sources (Wilson loop) to the corresponding Green's function in pNRQCD 



241 ]. The Green's function in pNRQCD describes the propagation of a static QQQ state in 



the colour representation C through a potential Vq- Loop corrections due to gluons of energy 
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FIG. 1: Static Wilson loop with edges x\ = (xx,Tw/2), x 2 = (x 2 ,Tw/2), 23 = (x 3 ,Tw/2), 
yi = (xi,—Tw/2), y 2 = (x2, — Tw/2), 2/3 = (#3,— Tw/2) and insertions of the tensors C^- k and 
C^ n in X = (R,T w /2) and Y = (R, -T w /2) respectively. 

and momentum of order a s /r contribute at next-to-next-to- next-to-leading order (N 3 LO) 
and are beyond the accuracy of this work. The matching condition valid up to and including 
NNLO is, for T w -> 00, 

(0\C u WC v ^\0) = Z c (t)exp(-iV c (t)T w )(0\S2 v (T w /2,-T w /2)\0). (1) 

The left-hand side stands for the expectation value of the three-quark static Wilson loop: 
a possible choice is shown in Fig. [TJ The static quarks are located in xi, x%, and 033 and 
propagate from the initial time — to the final time Tw/2. The colour tensors C u and 
C v ^ are inserted in the Wilson loop in the centre-of-mass coordinate R = (xi + X2 + Xz)/3 at 
the final and initial times respectively. |0) on the left-hand (right-hand) side stands for the 
vacuum state of QCD (pNRQCD). In the right-hand side, Vq stands for the static potential 
in the colour representation C, Zc for a normalization factor, Sc for the Wilson loop C U WC V ^ 
with all the quarks located in the centre of mass, and r = {ri, r 2 , r 3 } for the set of distances 
between the quarks, defined as 

ri = xi-x 2 , r 2 =xi-x 3 , r 3 = x 2 - x 3 ; (2) 

only two of these three distances are independent: I - ! + r 3 = r 2 . The explicit expressions of 
the three-quark static Wilson loop shown in Fig. [1] and of Sq v are 

C u WC vj < = C£ fc 0«/(R, xi, T w /2)<f>it r (T w /2, -T w /2, x 1 )0 rf (x 1 , R, -T w /2) 
x <p jj ,(R,x 2 ,T w /2)(j) j , s (T w /2, -Tvy/2,x 2 )0 sm (x 2 ,R, -T w /2) 
x ^(R,x 3 ,T w /2)^(T w /2,-T w /2,x 3 )^ n (x 3 ,R,-T w /2)C^ n , (3) 
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Sc v = C? jk cf>u(Tw/2, -T w /2, R)<p jm (T w /2, -TV/2, R)0 fen (T w /2, -T w /2, R)C^ n , (4) 

where repeated indices are implicitly summed from 1 to 3. The tensor C u is inserted at 
X = (R, Tw/2), while its conjugate C vjf is inserted at Y — (R, — Tw/2). The function <f> 
stands for a Wilson line: the spacelike Wilson line at time t reads 

0(y, x, t) = P exp ^ y ds (y - x) ■ A(x + (y - x)s, t)^ , (5) 

while the timelike Wilson line at position x reads 



(j)(t f ,ti,x) = P exp \ ig jf ' dtA°(t,x)j . (6) 

In both expressions, = A^T 0, and P stands for the path ordering of the matrices A^ along 
the Wilson line. In Eqs. ([3]) and (HI), we have explicitly written the colour indices of the 
Wilson lines in the fundamental representation. 

Let us now specify the colour representations C u . A QQQ state can be decomposed into 
the following representations: 

3<g>3®3 = 1©8©8©10, (7) 

where the singlet representation is totally antisymmetric, the decuplet is totally symmetric, 
and the two octets have mixed symmetries. A generic representation C u has three colour 
indices, i,j,k, running from 1 to 3 and is written in detail as C^ k . The labels u,v refer 
to the type of colour representations, specifically, when C and & are both in the singlet 
representation, the indices u and v are suppressed; when C u and C v ^ are both in the decuplet 
representation, u and v range from 1 to 10; when C u and C v ^ are in the antisymmetric or in 
the symmetric octet representations, the indices u and v range from 1 to 8. The concrete 
choice that we have operated for these rank-three tensors is given in Appendix El In the 
singlet and decuplet cases, C u and C"^ are real numbers. In the octet case, since the octets 
mix, it is more convenient to consider C u and as 2 component vectors; we will detail 
about this in the next section. 

The quantity (0\Sq v {T w /2, — T w /2)\0) is dimensionless. In perturbation theory, it may 
depend on Tw only logarithmically, therefore lim 1/Tw x hi(0\Sc v (Tw/2, — 7V/2)|0) = 0. 

Tw^oo 

Also lim 1/Tw x In Zc(x) = 0, because Zc(x) does not depend on Tw- Hence, the matching 

Tw— >oo 
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condition ([I]) may be rewritten as 

V ' T w ^oo T w C u Ctna 

where we have kept in the denominator a colour tensor normalization factor (cf. Eq. (1A5j) ). 
It is convenient to define 

<0| ^^r |0) = 1 + -M (0) (C, t) + t) + M®(C, r) + . . . , (9) 

with the quantities Ai^ n ' encoding all contributions of order g 2n+2 ~ for a given colour 
representation C. Analogously we may write 

V c (t) = V^°\x) + Vf\x) + vf (t) + . . . , (10) 

where V^ n \x) encodes all contributions of order g 2n+2 to the potential. From Eqs. (JSJ), (jUJ) 
and (fTUj) . the order by order matching conditions for the potential read 

M<°>(C,t), (11) 

fM^(C,t)-^M^ 2 (C,t)Y (12) 

>l (2) (C,r)->l (0) (C,r)>l {1) (C,r) + ^Al (0)3 (C,r)^ , (13) 

Note that the subtraction terms, M {0)2 ~ T w , M [0) M {1) ~ T w and .M (0)3 ~ T w , are diver- 
gent in the Tw — >■ oo limit. They cancel against divergences in Ai^ and .M^. Canceling 
the divergences may be interpreted as reconstructing the exponential exp (— iVc{x)Tw) in 
the matching condition ([1]). For this reason, the procedure of verifying the finiteness of the 
limits (|T2l) . (11311 . ... is often referred to as verifying the potential exponentiation. 





W 
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TV— >oo 
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III. THE STATIC POTENTIAL AT LO 



To set up the notation and to discuss the octet mixing, we start by calculating the three- 
quark static potential at LO, i.e. Vq . The calculation can be split into two steps: the 
computation of the amplitudes and the calculation of the colour factors, which will differ for 
each potential. Throughout the paper we choose the Coulomb gauge for the calculation of 
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FIG. 2: The three terms contributing to the static potential at order g 2 . Dashed lines are longitu- 
dinal gluons. 

the amplitudes, since it consistently reduces the number of diagrams to be computed. Of 
course, the calculated LO, NLO and NNLO potentials are gauge invariant 1 . 

At order g 2 , the diagrams that contribute to the potential are those shown in Fig. [2J 
in which a gluon is exchanged between two quark lines, thus leaving the third quark line 
untouched. In the following, we will call such a line a spectator line. Since diagrams involving 
gluon exchanges between quark lines and strings contribute only to the normalization Zc, 2 
we will adopt, in the following, a simpler representation of the diagrams without end-point 

3 

strings (see e.g. Fig. [3]). It is convenient to define MS°\C,t) = Vj.M^(C, r q ), where 

9=1 



1 Possible complications arising in Coulomb gauge because of the so-called Schwinger-Christ-Lee terms 
- 4^|. which involve a 2 -suppressed non-local interactions with transverse gluons, affect the potential at 
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next-to-next-to-next-to leading order or smaller and are beyond the accuracy of the present work. 
2 Our work is concerned with the static potential up to NNLO. Up to this order, two-body diagrams 
involving gluon exchanges with strings are of the same type as those encountered in the evaluation of 
the quark-antiquark potential and do not contribute to the potential by the same arguments used there 
B 44 1 (a more detailed discussion can be found in 45|). At NNLO order, there is also a class of 



three-body diagrams involving gluon exchanges with the strings. This class of diagrams has a transverse 
gluon emitted from one string and three longitudinal gluons coupled to it and to three different quark 
lines. These diagrams vanish because either they involve triple-gluon vertices with two transverse and 
one longitudinal gluon, but zero inflowing energy (taking Tw — > oo is equivalent to set to zero the energy 
flowing from the quark lines, see footnote [3]) , or they involve quartic-gluon vertices with one transverse 
and three longitudinal gluons. 
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A4 q (C,r q ) is the amplitude of the one-gluon exchange between two of the three sources: 
q = 1 corresponds to the exchange between the quark in xi and the one in x 2 , q = 2 
corresponds to the exchange between the quark in x x and the one in x 3 , and q = 3 corresponds 
to the exchange between the quark in x 2 and the one in X3. 3 The potential at order g 2 then 
reads 

^c (0) W= T lim fEM 0, (^,) = t/r(C)^ (14) 
Tw ^ool w ^ ^ |r,| 

The colour part of the amplitude has been factored in the colour coefficient fq (C). This 
coefficient is defined as 



nu -7-9(0) pv\ 

fWfQ) - jkl 'H'kk'W^j'k'V ^ 



where Tj$ww ~ T jj' T kk'^i' > Tffkk'W ~ T h^> kk ' T w and ^ffkk'W ~ ^ji' T kk' T w- 
In the singlet case, C = C* = S (see Eq. flAll) ) . we have 

fi°\S) = ti 0) (S) = f! > \S) = ~ (16) 

and the LO colour- singlet static potential has the well-known form: 

^ (0) W = -^fo + o + oV ^ 

3 Vl r il l r 2| N/ 

We note that in the limit where one quark is put at infinite distance the above potential 
should reproduce one of the two quark-quark potentials, either the antisymmetric antitriplet 
one or the symmetric sextet one. Since the singlet is antisymmetric one recovers indeed 



the antisymmetric triplet quark-quark potential 46] . Moreover, we observe that the colour- 
singlet quark- ant iquark potential is twice each quark-quark component of (fT7|) : we will 
generalize this result in Sec. [V] 

In the decuplet case, C u = C u * = A u (see Eq. (1 A4j) ) . we have 

/ r(A) = / f(A) = / f(A) = i (18) 

and the LO colour- decuplet static potential reads: 

^Mr + r + r> < 19 > 



3 For a baryon of N quarks in SU(iV), there will be N(N — l)/2 possible gluon exchanges, with (N — 2) 
spectator quarks. 
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FIG. 3: A diagram contributing to the mixing of the two octets at LO. The three quark lines 
represent from above to below the quarks in Xi, x 2 and X3 respectively. 



We note that in the limit where one quark is put at infinite distance the above potential 



46] 



reproduces the symmetric sextet quark-quark potential 

In the octet case, the one-gluon exchange mixes the symmetric and the antisymmetric 
octets, i.e. there is a nonzero colour amplitude with an initial symmetric octet state and 
a final antisymmetric one and viceversa. It is, therefore, convenient to define a potential 
Vo(t), which is a 2 x 2 matrix, and a vector colour representation: C a = 0^ k - ' 

A possible choice for the symmetric and antisymmetric octet representations OfA and Of^. 
is in Eqs. (IA3j) and (IA2j) respectively: in this choice, both representations are symmetric 
or antisymmetric in the first two indices i and j. According to the definition fl3]), the third 
index, k, is associated to the quark in x 3 , therefore we expect that the diagrams responsible 
for the mixing are those involving gluons attached to the third quark line, like the one shown 
in Fig. El Indeed, by computing Aif\0, r q ) it follows that the three 2x2 matrices fq(0) 
are given by 



A (0) (O) 







3 

§ 



ff(0) 




V3 
4 

12 



f?\0) 



1 




12 


1 


V3 


5 


4 


12 



(20) 



where we see that fi 0) (O) is diagonal but / 2 (0) (O) and / 3 (0) (O) are not. Hence, the LO 
colour-octet static potential reads 



V#>(t) 



1 

kll 



+ 







1 








V3 


12 








4 


^3 






5 


4 


St 




12 



I r 2 1 \ _xi L 



(21) 



The part of the potential proportional to l/|i"i| is diagonal and its entries are equal to the 
1/1 i"i I parts of the colour-singlet and colour- decuplet potentials. This can be explained by 
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observing that if the quark in X3 is put to infinity the two octets disentangle and we are left 
with two, antisymmetric and symmetric, quark-quark potentials. The parts of the potential 
proportional to 1 / |i"2 1 and 1 / J r 3 1 have the same diagonal elements but opposite off-diagonal 
ones: this means that they share the same eigenvalues but have different eigenvectors. The 
eigenvalues are (—2/3, 1/3) with corresponding eigenvectors 



A_ 2/3 = , O s ^j , A 1/3 = (±V30 A , O' 



(22) 



where the upper sign refers to the matrix with positive off-diagonal elements and the lower 
sign to the other one. If we construct a matrix P such that PMP~ X is diagonal, where 
M is one of the two non-diagonal matrices, the other being M', then neither PM' P~ l nor 
Pfi (0)P~ 1 are diagonal. The diagonalization of the part of the potential proportional to 
the distance l/|r 2 | (l/|r 3 |) thus simply corresponds to changing to a new octet representation 
symmetric and antisymmetric in the indices i and k (j and k). Note that by pulling at infinite 
distance the quark in x 2 or in xi we are left with a matrix, which, after diagonalization, 
reproduces again the two, antisymmetric and symmetric, quark-quark potentials. 

The fact that the two octets mix has, to our knowledge, not been discussed in the literature 



so far. The octet QQQ potential can be extracted from the lattice data in [40j . There, 
equilateral geometries (|ri| = |r 2 | = 1 1*3 1 ) have been taken into account for which the off- 
diagonal elements cancel (see Eq. ( 12T1) ). 4 In general, off-diagonal elements cancel in any 
isosceles geometry. Clearly, the mixing needs instead to be properly accounted for in any 
lattice simulation based on non-isosceles geometries. 



IV. THE STATIC POTENTIAL AT NLO 

The NLO, i.e. the order g 4 , contribution to the QQQ potential in the different colour 
representations is what we have called V^p. Two classes of diagrams contribute: two-body 
diagrams and three-body diagrams. These are shown in Fig. HI 

Two-body diagrams are simply the quark- ant iquark diagrams of order g 4 , which we know 
from the static quark-antiquark potential, with the static antiquark propagator replaced by 

4 For equilateral geometries, the singlet and octet potentials are attractive while the decuplet one is repulsive. 
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FIG. 4: Diagrams appearing at order g 4 in the three-quark potential. 

a quark propagator and with the addition of a spectator line. Their colour factor is of course 
different but the amplitude can be easily obtained from the QQ equivalent. 

Three-body diagrams such as the ones in Fig. H] d) and H] e) do not contain a spectator 
quark. We will show that diagrams of type H] d) only contribute to the exponentiation of 
the LO potential, i.e. cancel in Eq. f fl2|) against — A4^ 2 (C, t)/2, whereas the ones of type 
H]e), which include also diagrams with two gluons attached to the same quark line, vanish 
because they involve triple-gluon vertices of only longitudinal gluons. 

A. Calculation of V c {1) 

We start by examining the two-body diagrams in Coulomb gauge. These are shown in 
Fig. H]a), b) (the ladder and crossed diagrams), c) (the Abelian vertex correction) and f) 
(the gluon self-energy diagrams). In Coulomb gauge, the crossed diagram and the Abelian 
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vertex vanish. For instance, in position space the crossed diagram is proportional to 

7V/2 pTw/2 pTw/2 fTw/2 

dx° / dx'° / dy° / dy'° 9(x'° - x°)9(y° - y ,0 )5(x° - y°)5(x' - y'°) 

■Tw /2 J — /2 J —Tw /2 J ~T\y /2 
^T w /2 rT w /2 

dx° / d 

where the thetas come from the static quark propagators and the deltas from the longitudinal 
gluon propagators in Coulomb gauge. A similar argument applies to the Abelian vertex. In 
the case of the ladder diagram, the product of deltas and thetas does not yield zero but 
the complete result is 



/T w /2 fT w /2 
dx° / dx'° 6(x'° - x°)6(x° - x'°) = 0, 
■Tw /2 J - Tw /2 



2 



c? 3 q 



T 2 1 



^(Clw-f^. (23) 



(2tt) 3 q 2 J sJq y nM 2 r; 

where f q (C)i a d is defined as (we chose k and k! to label the spectator line) 

nu rparpb rpa rpb pv\ r 

i\ = ru n v\ • I 24 ) 

L 'mno L ' r nno 

In total, there are three ladder diagram contributions, ^4^(0, r g )i ad , with q running from 1 
to 3. 

Let's consider now the three-body diagram in Fig. H]d). We call M. q l q ](C, r q , ry) 3 body its 
contribution to Ai^(C, r), which is given by 

rT w /2 rT w /2 rT w /2 t>T w /2 

M§(C,T q ,T q ,) 3hody = -g 4 f^(C) 3hody / dx° dx'° dy° dy' 6(y°-x ) 

^ — Tw /2 J —Tw /2 J — Tw /2 J —Tw / 2 



x5(x'° -x°)S(y' -y°) 



d 3 q e iq r " f d 3 q' e iq ''V 



(2tt) 3 q 2 J (2tt) 3 q /2 
-^4lhCUo dy ^-^—, (25) 



In the case of Fig. H]d), the colour factor /^(C^body is defined as 



nu rp a rpb rpa rpb pVJ 

J qq/ (C) 3body t • (20 J 

L mno Lm ™ 



In total, there are six three-body diagram contributions, M qq ](C,r q ,r q >) 3hody , with q and q' 
(q 7^ <?') running from 1 to 3. 
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The contributions of the ladder and three-body diagrams cancel in Eq. (1121) against 



-A4 (0)2 (C,t)/2 



^MfiC^Yg) /2. This happens because 



Vg, 
Vg, 4- 



(27) 
(28) 



We will prove these identities for all representations C in the following Sec. IIVBI Hence 
the ladder and the three-body diagrams only contribute to the exponentiation of the LO 
potential. 

Finally, we are left with the evaluation of the diagram in Fig. H]f). The diagram has, 
in general, a q° dependence, however, the integration over time in the TV — > 00 limit sets 
q° = 0. 5 The fermionic part is gauge invariant; the gauge part, in Coulomb gauge, may 



be read, for instance, from 



471 | . The one-loop gluon self-energy contribution to the gluon 



propagator in momentum space and at q° = is 



i5 ab a. 



i 2 Air 



31 10 



2 \, AT 

11 -rn hi ^ + Y~Y nf , 



(29) 



where n/ is the number of massless light quarks contributing to the fermionic part of Fig. H]f ). 
The divergence has been renormalized in the MS scheme and /i is the renormalization scale. 
The contribution to the potential is 

3 



v c (1) w = E4 0) (e) 



q=l 

3 



d 3 q 



Ana* a s 



q 2 An 



11 



-n 



. n 31 

In —r -\ 

q 2 3 



10 



-n 



A-K\r q 



(30) 



where je is the Euler-Mascheroni constant, a\ = 31/3 — 10n//9 and {3q = 11 — 2rif/3. 

Since, in Coulomb gauge, all other diagrams of Fig. H] either vanish or contribute to the 
potential exponentiation, the contribution coming from the diagram in Fig. H] f ) is the only 
contribution to the potential at NLO. It has the same colour factor as the LO one, which 



5 The Tw — > oo limit comes from the matching condition ([5]). It sets q° 
dt / dq° exp(-iq°t) g{q°) = g(0). 

-T w /2 J 



as in 
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factorizes in front of the complete expression of the potential up to NLO. This reads 

^w = E4 0) (c)- 

q=l 



«s(l/K|) 



Qt 

1 + ^ (2/3 7£ + < 



(31) 



where the colour coefficients fq°\C) may be read from Eqs. (TT6l) . (TT8l) and (1201) . We recall 
that, in the octet case, Vo is a 2x2 matrix. 

The main outcome of Eq. ( }3~Tj) is that at NLO the QQQ static potential and the QQ static 

potential 44| just differ by the overall colour representation, but that the effective coupling 

r ol i — 

of the potential, ay(l/|rJ) = a s (l/|rJ) 1 + — - (2/3 7£ + ai) , is the same for all QQ, QQ 

and QQQ colour representations. There is no reason to believe that this result keeps holding 

at NNLO. Indeed, it has been shown in 48( that the colour-singlet and colour-octet effective 

couplings for the QQ potential differ at NNLO. 

In Feynman gauge, besides the diagram in Fig. H]f), also the diagrams in Fig. H] a), b) 

and c) contribute to the potential. The situation is very similar to the quark-antiquark case 

and it is straightforward to check that the final result up to NLO agrees with Eq. f lBTT) . 



B. Colour factors in the one- loop exponentiation 

In this section, we prove Eqs. ( 1271) and ( 1281) for all colour representations. In the singlet 
case, for all q and q' we obtain 

/ g (1) (5)lad = /^(5)3bod y = i (32) 

Together with Eq. ([TBI , this proves Eqs. ff2~T|) and ff28|) . Analogously, in the decuplet case, 
for all q and q' we obtain 



/i 1) (A)lad = /^(A)3bod y = ^, (33) 



which again, together with Eq. ( IT8|) . proves Eqs. (1271) and (|28|) . In the octet case, fq (0)i a d, 
as defined in Eq. f[2"lD . f^! (C)^ d y , as defined in Eq. (12"6T) . and fq°\0), as defined in Eq. 



201) . are 2x2 matrices. By explicit computation, one can show that 

/! (1) (0) la d= f 1 ° ] =(A (0) (o)) 2 , / 2 (1) (0) lad = ( * t 3 \=(ti°\o)f 





/ 3 (1) (0)iad = » *f \=(f?\0)T, (34) 
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which proves Eq. (J2 



One can also show that 



f (i) 

J 13 



(O) 



f(l) 

J 12 



3body 



3body 



f(l) 
i32 



3body 



This is enough to prove Eq. 




f?\0)f?>(0) 



(0), 



= /i (0) (O)/ 2 w (O), 

/f (O)/ 2 (0) (O). 
because, fq (O) being symmetric, 



(o), 



that(/< 0) (O)/®(O)) T = 4, 0) (O)/f(O), and, moreover, (f$(0) 3hody ) 



(35) 



/g'g(^)3body 



it holds 

6 



V. THE CASE OF N COLOURS AND N QUARKS 

A generalization of Eq. (|7|) to N > 3 colours but with three quarks can be easily obtained 

using the Hook length formula on the corresponding Young tableaux 49( . However, a system 

made of three quarks and N ^ 3 colours does not contain a colour-singlet state. For this 

reason, in the following, we will consider the case of N quarks and iV colours. With the 

increase in the number of quarks, also the number of representations increases rapidly, but we 

will always have a totally antisymmetric representation (the colour-singlet one) and a totally 

(2N - 1\ 

symmetric representation, whose dimension is , i.e. the number of independent 



N 

entries in a totally symmetric tensor of rank iV with indices running from 1 to N. 

The singlet representation ( 1A1I) can be easily generalized to any given number N of colours 
and quarks using the Levi-Civita tensor eyfc... of rank N. Since 

£ijk..Aijk... = N\ , (36) 

where repeated indices are summed from 1 to N, the normalized totally antisymmetric singlet 
tensor is given by 

~ ~\7m' 



The fact that the f^l (O)3t>ody matrices are not symmetric under the exchange q f-> q' does not contra- 
dict time-reversal invariance, since, in the complete three-body amplitude, for each diagram proportional 
to 4 } (0) 3body 

/(|rj |r g /|) there is a diagram proportional to |(0)3body/(|r g | |r g '|) that restores the 

symmetry. 
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where, from now on, a tilde will designate representations with iV colours and quarks. We 
provide now an expression for the colour factors fq°\S) relevant at LO. Since the singlet 
tensor is totally antisymmetric, the factors fq°\S) are equal for all q. The product of 
two Levi-Civita tensors can be expressed as a determinant of Kronecker 5 symbols in the 
following way: 



&im $in &il • • • 



^■ijk...^mnl. 



det 



3km Ofcn Okl 



\ 



(3? 



/ 



which generalizes the three-dimensional identity e ijk ei mn = S a (Sj m 8 kn -Sj n S km ) z^imi^jihn 



8j n 8ki) +8i n {5ji5km — 8jm$ki)- Using this property we obtain the colour factor [39] 



^■ijkl... 



rpa rpa c c 
1 im 1 jn°koOlp ■ ■ ■ 



'mnop.. 



N 



(39) 



where C F = (N 2 - l)/2/N. 

The singlet LO potential is then 



V- 



C h 



N 



n 1 
q=l yl 



(40) 



where the sum runs over all n = N(N — l)/2 possible one-gluon exchanges between two 
different quark lines and r is the dimension iV(iV — l)/2 vector (xi — X2, ...,xi — xat,X2 — 
X3, ...,Xjv-i — Xjv). We observe that the singlet QQ potential is — 1 times each two- 
body component of the singlet potential of a baryon made of quarks, which generalizes 
the well-known result that the quark-quark potential in an ordinary baryon (N = 3) is 
half the quark-antiquark potential. This may be understood in the following way: if we 
collapse N — 1 quarks in the same position the remaining one will "see" A^ — 1 times the 
quark-quark potential. This, in turn, corresponds to the quark-antiquark potential, since 
the SU(A^) antisymmetric representation of rank A^ — 1 describing a system of N — 1 quarks 
in a totally antisymmetric colour state has dimension A^ and corresponds to the conjugate 
of the fundamental representation, i.e. the representation describing an antiquark. 

For what concerns the totally symmetric representation, let A^ k be a generic symmetric 
tensor, with u running from 1 to ( j . The totally symmetric equivalent of Eq. (1551) 
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IS 

(2N-1 

N 



K jk ..Mf«... = M E ■ ■ ■ , (4i) 

n=l ' a(i>j'k'...) 

where, on the right-hand side, there are N Kronecker deltas and the sum is understood to be 
performed over all permutations of the indices i'j'k' .... The tensors A^ fc are normalized 

as 

^kL.Ahi- = 5UV ■ (42) 
In analogy with Eq. (|39|) . the totally symmetric colour factor relevant at LO is 

/\U rp a rpa XX AV fl 

f {0) (K) = jM " — J j ' _ • • • ^i'i'k'V... _ w / 4 on 

g " A\, AV.,, ~ N + I 1 ; 

ijkl... ijkl... 



The result follows from 



* 2N 2 ' 1 J 



the totally symmetric nature of Ay W and the normalization (jl2i) . The LO totally symmetric 
potential is then 

where, as before, n = ]V(jV — l)/2. 

We prove now the exponentiation of the colour-singlet and colour-symmetric potentials 
at NLO, i.e. Eqs. ( 1271) and ( 1281) . for a baryon in SU(iV) made of N quarks. We can write 
the colour factor fq (S)^ as 

fWfQ} _ e ijkl... rpa rpb rpar-pb c c E-mnop. 
in Wjlad /— -rr 1 ix 1 xm jy yn°ko°lp ■ • • 



^ ix 1 xm 1 jy 1 yn°koVlp ■ ■ ■ ' 



and the colour factor / g (A)^ as 



An rparpb rpa n^b XX All 

Jq l^Jlad 



A".,, AV.,, 

tjkl... tjkl... 

Using Eq. (PHI) , the totally antisymmetric nature of e mnop ,,,, the totally symmetric nature of 
AV. fcL and the normalizations (|36|) . (j42l) we obtain 

2 

.iVTl 1 



fM(S,AU= (-§-) V,/. 
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where the upper sign refers to the antisymmetric case and the lower sign to the symmet- 
ric one. This proves that fg (S)^ and fg (A)i ad are the squares of fq (S) and fq°\A) 
respectively, i.e. Eq. ( j27l) . 

For the three-body diagram we adopt a similar procedure, with the difference that here 
the contracted indices will be iV — 3. 7 The colour factors are then, for all q and q ; , 



and 



fqj {S)3body 



fql> ( A) Sbo&y 



rpa rpb rpa rpb £ 

1 ix 1 xm 1 jn 1 ko°lp 



A« Ay 

ijkl... ijkl.. 



Proceeding like before, we obtain 



fqt>(S> A)3body 



iV =f 1 



'mnop.. 



<N\ 



Au i~pa rpb i~pa rpb X /\v 
^ijkl.. 1 ix 1 xm 1 jn 1 ko°lp • • • ^rr 



mnop. . 



Vg, q', 



(46) 



which proves Eq. 
case. 



for the antisymmetric (upper sign) and the symmetric (lower sign) 



VI. THE THREE-BODY PART OF THE STATIC POTENTIAL AT NNLO 

We may ask when a genuine three-body interaction, i.e. a contribution which is not 
the sum of three l/|r 9 | terms and is not generated by the exponentiation of two-quark 
interactions, shows up in the Wilson loop. This happens at order g 6 . More precisely, we 
write 8 

3 2body//n 

r ( 2 ) (r) = Vc 3body (r) + a 3^^0 ) (4?) 

q=l \ T «\ 

where the three-body part of V c f , V^ 3body , is defined as the part of V^ 1 that vanishes when 
putting one of the quarks at infinite distance from the other two, i.e. in the limit |rj|, 



7 For definiteness, we assume the two gluons to be attached to the same quark line. However, starting from 
N = 4 quark lines, it is also possible that a gluon is exchanged between two quarks and a second one is 
exchanged between two different quarks. This is again a l/(|r 9 ||r 9 '|) term and by similar arguments it can 
be shown that its colour factor is also the square of (|5TJ|). thus obeying 



8 We assume that ln(/i|r 9 |) terms have been resummed such that the potential up to NLO reads 



3 



Ve(t) = >J U) (C) 



q=l 



L <-ll 



Q! S (1/K1) . 
1 H -. (2/?o7b + ai ) 

47T 



Under this condition, terms like ln(/x|r g |) 



or ln 2 (/i|r 9 |) are absent at NNLO. 
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FIG. 5: Three-body diagrams at NNLO that contribute, in Coulomb gauge, to the exponentiation 
of the LO and NLO potential. 

\Tj\ — > oo [i 7^ j) with fixed \rk\ {k ^ i and k ^ j). Since Vq is gauge invariant, then, by 
definition, also the numerical coefficients a^ hody (C) and V^ 3body are. V^ 3body may only stem 
from diagrams with gluons attached to all three quark lines. 

At order g 6 , we have many diagrams that involve gluons attached to all three quark lines. 
These can be divided into some basic categories. The adoption of the Coulomb gauge proves 
again useful, making only a small subset of these diagrams different from zero. We thus 
have the following diagrams, evaluated, for simplicity, between totally antisymmetric and 
symmetric colour states only. 

1. The diagrams displayed in Fig. [5] contribute to the exponentiation of the tree-level and 
one-loop potentials. At this order of perturbation theory, the matching condition is 
given by Eq. ( fT3l) . It is easily shown that the amplitudes of the diagrams a), b) and 
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c) are 9 



T 3 1 -* r T 3 

-m — rrn> — r*~ 
-«l l r g'l l r <?l — 



o\ r„ r„' r„ — o: r„ r„/ 



^ = ^/; y (4 i r |i r !iir,r (49) 



Keeping in mind that there are three diagrams of the form of Eq. (|48p , one of type a) 
and two of type b), for each q, q' pair, and six diagrams of the type of Eq. (j4"9"|) . it is easy 
to see that their contributions cancel against -A4 (0) (C, x)M w (C, r) + A4 (0)3 (C, t)/3 
in the matching condition flT3|) and therefore do not contribute to Vq 2 \x). 

The amplitude of diagram d) can be obtained from Eq. (|25|) substituting one of the 
two longitudinal gluon propagators with Eq. ( )29l) . yielding 




= [2A,(ln(Ai|r fl |) + 7*) + «i] 4 1 ' ) (C)3bo dy ^frr^n- (50) 

I 9 1 I <l' I 

Recalling that we already proved the exponentiation relation of the colour factor 
/gi (C)abody i n Sec. IIVBI from Eqs. (jHj) and (13T1) we see that diagrams of the type d) 
(two for each q,q' pair) cancel against —M^(C,x)M^\C,x) in Eq. ( TT3l) . 

2. Abelian diagrams such as the ones in Fig. [6] are easily shown to be zero in Coulomb 
gauge. However in different gauges, such as the Feynman gauge, these diagrams are 
expected to give a contribution to the exponentiation and a contribution to the order 
a 3 result, as their two-body counterparts do in the QQ case js, 9, 44]. 



3. The non- Abelian diagrams shown in Fig. [7] also vanish. Diagram a) has a vanishing 
colour factor between singlet- singlet and decuplet-decuplet initial-final states and, in 
Coulomb gauge, a vanishing amplitude as well. The dashed blob in diagram b) is a loop 
of gluons or fermions. Lorentz invariance dictates that its Lorentz tensor structure has 



In the general case of a baryon in SU(iV) made of N quarks, 



A) = fggi(S, A) = f c qqlq „(S, A) = ( t 



C F n ' 



for all q, q' and q" . The upper signs refer to the antisymmetric (singlet) representation and the lower signs 
to the symmetric representation. 



22 



d) 



e) 



f) 




FIG. 6: Abelian three-body diagrams that have zero amplitude in Coulomb gauge. 

to be composed by combinations of a metric tensor and the external momenta q^. 
Since the sources are static, this guarantees that the Lorentz structure is proportional 
to at least one power of qf. By means of the usual argument, in the TV — > 00 limit, q® 
gets multiplied by S(qf) and vanishes. Finally, also diagrams c) and d) vanish because 
they involve non-Abelian vertices with longitudinal gluons only. 

4. We are then left with diagrams of the type shown in Fig. |HJ in Coulomb gauge, these 
are the only ones contributing to the three-body interaction. 

We now proceed to the evaluation of the diagrams of Fig. [HJ There are six different 
diagrams of each type a) and b): for each source line there are two diagrams where this line 
couples to two gluons (like the bottom line in Fig|8]). These two diagrams are symmetric with 
respect to a permutation of the other two lines. We call T-Lq and T-Lq the momentum-space 
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a) 



b) 



c) 



d) 



FIG. 7: Non-Abelian three-body diagrams that have zero amplitude in Coulomb gauge. 



b) 



FIG. 8: The only three-body diagrams that are not exponentiations and that have a non-vanishing 
amplitude in Coulomb gauge. Dashed lines are longitudinal gluons, curly lines are transverse ones. 

amplitudes of the diagrams in Fig. [8] a) and b) respectively and "He = "He + ^c- 

We consider now the colour structure of the diagrams in Fig. [HJ The colour factors fn(S) 
and fn(A) are equal for all twelve diagrams: 



US) = - 1 - and /„(A) = ~ 



(51) 



We note that the singlet and decuplet colour factors share the same sign, hence also the 
contributions to the potential from these diagrams will share the same sign, at variance with 
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the tree-level and one- loop results. 10 

We compute now %c- We call q 2 and q3 the momenta that flow out of the first and 
second quark line. Setting to zero the external energies, we obtain 

f H (C)g 6 f d*k 4(q 2 -kq3-k-q 2 -q 3 ) 



and 



^c(q 2 ,q 3 ) 



^(q 2 ,q 3 ) 



q 2 q 2 J (27r) 4 (A;° + 2 e )(k-q 2 ) 2 (k + q 3 ) 2 (£; 2 + 2 e ) 
fn{C)g 6 f d 4 k 4(q 2 -kq3-k-q 2 -q 3 ) 



q 2 q 2 J (27r) 4 (-A; + 2e)(k-q 2 ) 2 (k + q3) 2 (P + ^) 



(52) 



(53) 



Summing and 1-L h c yields 



"Hc(q2,qs) 



fn(C)g 6 f d 3 k 4(q 2 -kq3-k-q 2 -q 3 ) 



qloi 



8q|q| 



(2tt)3 k 2 (k-q 2 ) 2 (k + q 3 ) 



+ 



IO2IIO3I 



q 2 • 


q 3 + 


02 1 


lOal 


I02I 


loal 


|q2 + qsl 



I02I I03I 



The contribution of this diagram to the potential in position space is 



^ic( r 2, r 3 ) 



d 3 q 2 



3 «q2-r2 



d 3 q 3 



e iq3 ' r ^c(q 2 ,q 3 ); 



(2tt) 3 J (2tt) 3 

the total contribution of all twelve diagrams of the type shown in Fig. [8] is 

V$(r) = 2 [V nc (r 2 , r 3 ) + V wc ( ri , -r 3 ) + V nc (-r 2 , -n)] . 



(54) 



(55) 



(56) 



As shown in App. |Bj V^(r) may be expressed as a double integral suitable for numerical 
evaluation. We have considered the following geometries. 



10 For the antisymmetric and symmetric representations of a SU(iV) baryon made of N quarks, the colour 
factors are given by 

f„,(q\ E ijkl... rpd r 



r I q\ E ijkl... rpd rpa rpb rpe fbdcraecc ^mnop 

JH\d) — — 7^T 1 im, 1 m 1 kr I roJ J °lp ■ ■ ■ 



N\ 



MA) = 



A u rpd rpa rpb rpe rbdc-faecz Au 
LA ijkl... J -im 1 jn ± kr 1 roJ J °lp ■ ■ • *-*rr 



ranop. . 



A u Av 
ijkl... mnop.. 



Using (without summing over u), 

Eijkli ...In-3 ^mnoii---(jV-3 $im {3jn&ko $jo$nk) $in {Sjm^ko $jo$km) $io {.&jn$km $jm$kn) 



A U \U 

ijkli . . mnoli . . ./ jv-3 



we obtain 



N{N - 2) (AT- 1) 



f H (S,A) = - 



N{N + 2){N + 1) 
iV±l 
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(A) Isosceles geometry in a plane 
In this geometry, the three quarks are placed in different positions of the same plane, with 
two distances chosen to be equal: 1 1~2 1 = l r 3| = r and ?2 • £3 = cosd. The quarks are located 
at the vertices of an isosceles triangle. The potential V^q depends on r and 9; it has the 
form 

Vrc(r,6) = MC)af-^. (57) 



In Fig. 9(a), we plot c-n(9) as a function of 9. The coefficient is always positive, giving rise 
to an attractive contribution to the potential, both in the singlet and decuplet channels (we 
recall that the colour factors (l5Tj) are negative). The dependence on the angle 9, i.e. on the 
geometry of the configuration at fixed r, is weak: c%{9) ranges from a maximum of about 
1.46 at 9 0.65 to a minimum of about 0.49 at 9 = 71. On the contrary, the dependence 
on the geometry of the two-body contributions to the potential, such as Eq. f fl4l) . is much 
stronger. In particular, the two-body contribution diverges in |ri| = 0, i.e. for 9 = 0. 

The weaker dependence on the geometry of the three-body contribution with respect to 
the two-body contribution could signal the onset of a smooth transition towards the long- 
distance Y-shaped three-body potential seen in the lattice data. This long-distance potential 
turns out to depend only on one length, L, which is the sum of the distances between the 
Fermat point of the triangle made of the three quarks and the three quarks. For isosceles 
triangles, L has the following dependence on r and 9: 

cos - + v^sin § for < 9 < 2vr/3 
2 for 2tt/3 < 9 < n 

Note that the Fermat point of any triangle with an angle greater or equal than 2tt/3 is 
located at the vertex of that angle. In terms of L, Eq. fl57|) becomes 

Vr c (L,9) = MC)af-^l. (59) 



L = g(9)r, where g{9) = { 2 2 ~ ~ (58) 



In Fig. 9(b), for completeness, we plot g(9)c-u(9) as a function of 9. The plot is qualitatively 
very similar to the plot of c-h(9): the maximum gets shifted to 9 w 1.047, numerically 
equivalent to the equilateral geometry 9 = ir/3, which thus appears to be the energetically 
favored one for V^(L, 9) at fixed L. 

(A.l) 9 = 0: two quarks in the same position 
A special case of isosceles geometry is 9 = 0, where two quarks are located in the same 
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(a) 



(b) 



FIG. 9: In Fig. (a), we plot the coefficient c%{0) as denned in Eq. ([57]) and obtained from the 
numerical integration of Eq. (]B3p . In Fig. (b), we plot g(9)cu(9). 



position. From J d(q 2 • q 3 ) Hc(q2,q3) = 0, it follows that V H c(0, r 3 ) = V^c(i*2, 0) = 0, 
hence V^(r, 0) = 2Vuc( r , r )- The three-body potential is finite and given by: 



oh(0) 



with c%(0) = 6 — 



7T 



(60) 



(A. 2) 9 = 7r/3: planar equilateral geometry 
In the equilateral case, we have c^(n/3) ~ 1.377. We may compare the relative magnitude 
of the three-body contribution to the tree-level potential. In the singlet case (cf. Eq. f[T7]) ). 
the ratio yields 

V (°)(r) 4 a *W r >~ 2.90 ' 1 j 

where we have made explicit the scale dependence of the coupling constant. We note that, 
using a s at one loop, V^(r) may become as large as one sixth of the tree- level Coulomb 
potential in the region around 0.3 fm, where, at least in the QQ case, perturbation theory 
still holds Q- 

(B) Generic geometry 

In the most general geometry, the three-body potential ( 155]) depends on two coordinates. 
We may arbitrarily chose one of these coordinates to be L, leaving the other unspecified. If 
we call a, b, c the lengths of the three sides of the triangle made of the three quarks, then L 
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FIG. 10: The normalized three-body potential, V^(L, ...)/(— /%(C)ag), plotted as function of L 
for the geometry described in (B.l). 



is given by [37] 



a 2 + b 2 + c 2 a/3(o + b + c)(-a + b + c)(a - b + c)(a + b - c) 

Lj — ~r 

2 2 

2tt 

L = a + 6 + c — max(a, 6, c) for ^ > -r> 



for 9, 



max 




(62) 



where # max is the largest angle of the triangle. 

(B.l) Planar lattice geometry with two fixed quarks 
In Fig [TUl we plot the three-body potential obtained by placing the three quarks in a plane 
(x, y), fixing the position of the first quark in (0, 0), the second one in (1,0) and moving the 
third one in the lattice (0.5+0.125 n x , 0.125 n y ) with n x e {0, 1, ...,20} and n y e {0, 1, ...,24}. 
The plot clearly shows the dependence on the geometry at fixed L, however, the dependence 
is weaker than in the two-body case. 

(B.2) Three-dimensional lattice qeometry with the three quarks moving along the axes 



In the lattice calculation of Ref. [28], the three quarks were located along the axes of a 
three-dimensional lattice, namely at (n^O, 0), (0,^,0) and (0, 0, n z ), with n x 6 {0,1, ...,6} 
and n y ,n z G {1, 6}. For the sake of comparison, we consider the same geometry and plot 
the corresponding three-body potential in Fig. [TTJ The plot shows a weak dependence on 
the geometry: much weaker than in the two-body case, but also somewhat weaker than in 
the geometry considered in (B.l). 

As a final remark, we would like to note that V^, the contribution of the diagrams 
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FIG. 11: The normalized three-body potential, Vj£}i(L, ...)/(— /%(C)ag), plotted as function of L 
for the geometry described in (B.2). 

shown in Fig. E] calculated in Coulomb gauge, has an unambiguous physical meaning. From 
Eq. ( 1B3|) . it can be seen that this contribution vanishes when one of the quarks is put at 
infinite distance from the other two. Hence no two-body contribution gets entangled in V^, 
which can be rightfully identified with the three-body potential, V^ 3body , defined in Eq. (147|) . 

VII. THE COLOUR-SINGLET STATIC POTENTIAL AT NNLO 

In the colour singlet case, Eq. ( 14T1) becomes 



The coefficient a 2hody (S) is independent of the geometry of the three quarks: we can take 
advantage of this fact and calculate a 2hody (S) without performing any explicit two-loop 
calculation. In a configuration like the one described in (A.l), Vy is only a function of the 
distance r between one quark and the other two located at the same point: 



up to a singular term independent on r that we may drop, for instance, by dimensionally 
regularizing the potential in momentum space. In this configuration, V s (r) is equal to 
the static quark- ant iquark potential, because, when three quarks are in a colour-singlet 
configuration and two of them are located at the same point, these two behave as a an 




(63) 




(64) 
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antitriplet in colour space, i.e. as an antiquark. Owing to the two-loop result of the quark- 
antiquark potential, we may therefore write {9] 



Aai 1 



7T 



«2 + + ±1e) 06 + 1e (4ai/?o + 2/3i) 



3 r (4tt) 2 
where (3\ = 102 — 38n//3 and 

4343 2 9 , _ , /1229 52 > , A 100 a 

fl2 = it + " 4 + c( 3) ~ \w + y c(3) ) nf + ^r n '- 

From Eqs. (jM]) and fl65l) . it follows that 
2 1 



2body 



(S) 



3 (4?r) 



a 2 - 367T 2 + 3tt 4 + — + 4Ye) fit + 1e (4a x /3 + 2/?i 



(65) 



(66) 



(67) 



The complete NNLO expression of the three-quark colour-singlet static potential, Vy (r), is 
then given by Eq. f )63|) . where V s 3body (r) = V^°*(r) can be read from Eqs. fl56|) and (IB3|) . and 
a 2body (5') from Eq. ( 167|) . The explicit expression of the colour-singlet static potential up to 
NNLO is listed in Eq. ABB). 



VIII. CONCLUSIONS 



We have studied the static potential of a three-quark system in perturbation theory up to 
NNLO. Up to NLO, we have analyzed all the colour channels (singlet, octets and decuplet) 
of the SU(3) case and the results have been generalized to SU(iV) with iV quarks for the 
totally antisymmetric and totally symmetric channels. At LO, the potential is a sum of three 
Coulombic one-gluon exchanges between two of the three quarks. We have pointed out that, 
already at this order, octets mix. At NLO, after proving the potential exponentiation, the 
potential turns out to be simply a sum of two-body contributions, whose effective coupling 
ay is independent of the considered colour state and is the same as for the QQ, QQ and 
QQQ potentials. It is expected that ay becomes dependent on the colour state at NNLO, 
as it happens in the QQ case. 

At NNLO, the first genuine three-body contribution appears. Three-body contributions 
are specific features of the QQQ potential and for this reason of particular interest. We 
have calculated this contribution, providing numerical results for several geometrical confi- 
gurations. The general outcome is that the dependence on the geometry of the three-body 
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force is weaker than for the two-body force. Combining the three-body contribution with the 
two-body contribution extracted from the NNLO expression of the quark- ant iquark static po- 
tential, we have obtained the complete three-quark colour-singlet static potential at NNLO. 
It reads 



K(r) 



2^ a B (l/|r g |) f x | a„(l/|r g |) 
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(68) 



+16?rrVrV 



with R = xr 2 — yr 3 , R = |R| and 



, f 1 , R*R J 
x / ax / ay 

'0 io R_ 

M = \r 2 \^/x(l — x) + 1 1*3 1 a/ y(l — 1/). Note that by pulling one of the quarks at infinite 



distance from the others, the three-body potential as well as two of the two-body poten- 
tials vanish and Eq. fl68|) reduces to the quark-quark antitriplet static potential at NNLO, 
relevant for QQq baryons. 



In 



24j , also the three-loop leading logarithmic contribution in the infrared cut off has been 



calculated. Since that calculation does not account for the octet mixing, its result applies for 
geometries where the mixing cancels, like the isosceles one. It would be interesting to extend 
that calculation to generic geometries and combine the result with the complete NNLO result 
given above. 

Other possible future developments include comparisons with lattice results. They exist 
both for the ground state (the colour-singlet state) and for the possibly first gluonic excitation 
of the QQQ system 27J,l28|. An accurate comparison in the short range will show the running 
of the three-body potential and determine at which distances a perturbative description of 
the three-body potential breaks down. It may also serve to establish the nature of the 
gluonic excitation seen in the lattice data, determine if it is indeed the first excitation and 
clarify if, in the short range, the three static quarks assume a singlet, an octet or a decuplet 
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colour configuration; it may also serve to extract the masses of the gluelumps made of three 
static quarks. For all this it is crucial that octet mixing is properly taken into account 
in the analysis and in the lattice set up if geometries different from the isosceles one are 
used. Finally, in the case of more general geometries, it would provide particular insight 
in the non-perturbative dynamics of QCD, to investigate the transition region from (short) 
distances dominated by two-body forces (where the potentials depend on two coordinates) to 
(long) distances dominated by three-body forces (where, for the Y-shaped configuration, the 
potentials depend only on one string length). In this respect, the weak dependence on the 
geometry shown by our results for the leading perturbative three-body contribution could 
indicate a smooth transition to the Y shape. 

The QQQ static potential at higher order is relevant for the determination of the masses of 
the baryons made of three heavy quarks. Our NLO result is sufficient to provide the masses 
at NLO, 11 while at NNLO also 1/m and 1/m 2 potentials should be included. Clearly, having 
a reliable determination of the masses is of valuable help in the experimental searches. 



In 52J , the possible relevance of baryonic states in the quark-gluon-plasma phenomenology 
was pointed out and in [40] finite temperature lattice QCD simulations of QQQ systems in 
all colour channels were performed. The lattice data are very accurate also in the short range 
and clearly distinguish (in an equilateral geometry) among the singlet and octet (attractive) 
potentials and the (repulsive) decuplet potential before screening sets in. Temperature effects 
at short distances may be systematically included along the lines developed in Ref. {5^ for 
the QQ case and comparisons with finite temperature data may be performed. 

In general, one expects that QQQ states in a thermal bath will experience a much richer 
phenomenology than QQ states. First, more colour configurations are possible, second, 
among these, not only the singlet but also the octet states are subject, at least in some 
geometries, to an attractive interaction. Finally, there will be a larger variety of possible 
transitions among the different states induced by the thermal bath. Thermal transitions 
between colour- singlet and colour-octet or colour- decuplet states will likely be the dominant 



11 If implemented, our result may affect the mass determinations obtained in Ref. 50| within a variational 
study of weakly-coupled baryons. We note that the value obtained there for the bbb ground state is very 
close to the lattice determination of Ref. [51], providing an indirect evidence in support of the Coulombic 
nature of the system. 
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source of the QQQ colour-singlet thermal decay width in the short distance, low temperature 
regime as it is the case for the colour- singlet to colour-octet transitions in the QQ case [531.154 1. 
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Appendix A: Representations 



To ease the reader, we reproduce here from [24J the tensors for the singlet, two octet and 
decuplet representations in which the product of three triplet representations of SU(3) may 
be decomposed. The totally antisymmetric singlet tensor is 

$ijk — ~J=, (Al) 
the octet antisymmetric in the indices ij is 

D Aa — r-- T a r\Aa* _ t« ( K9\ 

^ijk ~ c ^]q ± kqi ^ijk ~ ^ijq- 1 qkt K^^J 

where the index q is summed from 1 to 3, the octet symmetric in ij is 



(~)Sa __ ' / rpa 

a/3 9 


+ e jkqT iq j , 


f)Sa* 
U ijk ~ 


1 

= 7I 


(tikqTqj + £jkqT qi ) 




and the symmetric decuplet is 
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One can easily check the following normalization and orthogonality relations: 

Q Q _ 1 f\Aa*f)Ab _ f\Sa*f\Sb _ cab act a p _ rap 
'Jijk'Jijk — -L, ^ijk ^ijk — W ijk W ijk — > ^ijk^ijk — > 

C D Ab Q Q A CT /^iAa*/~\Sb s^Aa* act ^\Sa* act n / ac\ 

>~>ijk^ijk — ^ijk^ijk — D ijk^ijk ~ w ijk w ijk ~ w ijk ^ijk ~ w ijk ^ijk ~ u - K 1 ^ ) 
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Appendix B: Position-space three-body potential 



From Eqs. (154"]) and (1531) . V^c( r 2> r 3) ma Y be written as 
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In order to evaluate the integrals, it is convenient to introduce the Feynman parameters x 
and y: 
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In this form, the integrals in q 2 , q 3 and k can be performed analytically and V^c( r 2, 1*3) ends 
up as a two-dimensional integral in x and y: 
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where R = xr 2 — yr 3 , R = |R| and M = |r 2 | \J x{l — x) + |r 3 | \/y(l — y). 
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